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Abstract 

We have used calculations based on density functional theory to investigate the energetics of hydrogen 
absorption in calcium-intercalated graphites. We focus particularly on the absorption energy and the 
stability of the hydrogenated material with respect to decomposition into graphite and calcium hy- 
dride, which is essential if this material is to be useful for practical H2 storage. The calculations are 
performed with two commonly used approximations for the exchange-correlation energies. Our calcu- 
lations confirm earlier predictions that the absorption energy is approximately —0.2 to —0.4 eV, which 
is favourable for practical use of Ca-intercalated graphite as a hydrogen storage medium. However, we 
find that the hydrogenated material is strongly unstable against decomposition. Our results therefore 
explain recent experiments which show that H2 does not remain stable in CaC6 but instead forms a 
hydride plus graphite. 
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1. Introduction 

Finding an alternative, or secondary, energy 
source to replace petroleum for use in vehicles is 
of paramount importance. Hydrogen is a possible 
candidate and, with the advancement of proton 
exchange membrane (PEM) fuel cell technologies 
|H, could provide a viable energy solution. The 
practical ambient storage of hydrogen has long 
been a challenge @, 0|. If hydrogen fuel cells 
are to be used in the transport industry then a 
lightweight energy dense storage medium needs 
to be found jij. The U.S. Department of Energy 
state that a storage medium needs to store 6.5% 
j^] hydrogen by weight at ambient temperature 
and pressure. In order that the material can take 
in and release hydrogen at ambient temperature 
and reasonable pressures, the absorption energy 
also needs to be in the range 0.2—0.8 eV (1.000 eV 



= 96.485 kJ/mol). 

Graphite intercalation compounds (GICs) have 
shown promise as a potential ambient high- 
density storage medium because they have a low 
density. They have the great advantage that their 
chemistry and the spacing between the carbon 
planes can be tuned p, 0] • They are also inexpen- 
sive to produce. Many GICs have been considered 
for hydrogen storage and possess very different 
properties when hydrogen is introduced into them 
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Recently calculations based on density 
functional theory (DFT) have been used to survey 
the binding energies of hydrogen in a wide range 



of metal-intercalated graphite materials The 
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results showed that alkaline-earth intercalates ap- 
pear to be very promising, because their calcu- 
lated binding energies are in the required range. 
An essential requirement for the practical use of 
any of these GIC materials is that they should be 
stable against decomposition into the metal hy- 
dride and graphite. Very recent experiments [12| 
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on the hydrogenation of calcium graphite (CaCe) 
suggest that this may be a major issue and that 
the attempted hydrogenation of the material leads 
to the formation of orthorhombic calcium hydride 
(CaH2) and the complete deintercalation of the 
starting CaC6 leaving graphite. This has also 
been shown to occur for lithium graphite (LiCg) 
(ljf . Surprisingly, this crucial question of stabil- 
ity against decomposition does not appear to have 
been addressed before by modelling. The recent 
DFT work by Cobian and Iniguez 11] on hydro- 
gen storage in GIC materials did explore the ques- 
tion of stability against metal clustering, but did 
not consider the possibility of hydride formation. 

The purpose of the present paper is to use 
electronic-structure methods to study the sta- 
bility of hydrogenated Ca-intercalated graphite 
against decomposition into graphite and CaH 2 . 



The experiments of Srinivas et al. |l2j indicate 
that this decomposition occurs spontaneously, so 
that it must be exothermic. This means that the 
Ca-GIC could not be used without modification 
as a hydrogen storage material. However, it is still 
possible that the hydrogenated material might be- 
come stable against decomposition if other atomic 
or molecular species were co-intercalated together 
with Ca. In considering whether and how this 
can be done, it is important to know the heat 
of decomposition of hydrogenated Ca-GIC. This 
important quantity is not available from the ex- 
periments, and our aim here is to use DFT to 
calculate it. 



In the next section, we summarise the essen- 
tial background information about Ca-GICs, and 
we also outline the DFT techniques used in this 
work. In section [31 we present our calculations, 
first on Ca-GIC (13. 2p itself at various stoichiome- 
trics, then on the hydrogenated material ( 13. 3p . 
and finally on the CaH 2 crystal (13.41) . These cal- 
culations yield values of the heat of decomposition 
of hydrogenated Ca-GIC, computed with different 
DFT approximations. As expected, the decompo- 
sition is exothermic, and we discuss in section H] 
the implications of our results. 



2. Background information 

2.1. The materials 

The structure of graphite is well known to con- 
sist of planes of carbon atoms, with each plane 
having a hexagonal honeycomb structure. Dif- 
ferent registrations of the carbon planes relative 
to each other are possible, the main two being 
A—B—A (hexagonal, 2H) and A—B—C (rhombo- 
hedral, 3R), with the 2H registration being most 
thermodynamically stable at standard tempera- 
ture and pressure The nearest neighbour 
C-C distance is 1.42 A. The layers are held to- 
gether by weak van der Waals dispersion, a phys- 
ical effect that is not included in standard forms 
of DFT. Nevertheless, it is known that the lo- 
cal density approximation (LDA) of DFT fortu- 
itously predicts an equilibrium interplanar spac- 
ing (3.38 A) that is very close to the experimental 
value (3.36 A). 

For Ca-intercalated graphite the registrations 
of the graphite layers are shifted so that car- 
bon atoms superimpose each other producing an 
A — A — A stacking. It has been shown by x- 
ray diffraction experiments [l5| that the Ca in- 
tercalants are positioned directly above and be- 
low hexagon centres (see figured]). The only ex- 
perimentally observed stoichiometry is C&Cq in 
a — (3 — 7 stacking (where we use the notation of 
Dresselhaus and Dresselhaus j^l with Greek let- 
ters labelling successive intercalant layers) such 
that each successive Ca layer is shifted relative 
to the preceding one resulting in a rhombohedral 
structure. This only exists as a stage I GIC, 
where the staging numbers how many graphite 
layers exist between each intercalant layer. It is 
known that other GICs produce different stage 
GICs, in the case of KC24 a stage II GIC. The in- 
terplanar spacing of Ca-GIC is increased with re- 
spect to graphite to 4.524 A (a 35% increase over 
graphite) and the C-C in-plane nearest neighbour 
distance to 1.444 A (a small 1.7% increase over 
graphite). 

A definable crystal structure for the hydro- 
genated Ca-GIC is not observed experimentally 



161 ] . Once hydrogenated the CaC 6 diffraction 



peaks weaken in intensity and the reflections from 



the CaH2 structure become apparent. A weak 
unidentified peak appears but is too weak to re- 
late to a major structure change. 

The CaH 2 crystal has an orthorhombic unit 
cell with Pnma symmetry (International Tables 
group 62). The unit cell consists of four CaH2 
units at each of the Wyckoff 4c positions. This 
means that the crystal requires only the x- and z- 
coordinates of one Ca atom and of two inequiva- 
lent H atoms to fully describe it; the y-coordinate 
is fixed by symmetry. The Wyckoff 4c positions 
are given in table [3j Experimental structural pa- 
rameters of the crystal have been reported by Wu 
et al. [17|, and calculations on the structure have 
been performed by Li et al. 18| and Wolverton 



and Ozolins 19 



2.2. Electronic- structure techniques 

This work employs standard electronic- 
structure methods based on DFT, applied to 
periodically repeated simulation cells. At the 
temperature of interest, electronic excitations 
are negligible, and all the predictions are derived 
from the electronic ground state energy of the 
system and the forces on the ions for given 
ionic positions. We use the projector augmented 
wave (PAW) method to determine the electronic 
ground state, and all the calculations are per- 
formed with the VASP code. We have always 
made efforts to ensure that the calculations 
are well converged with respect to the PAW 
plane-wave cutoff, so as to ensure basis-set 
completeness, and fc-point sampling to reduce 
or eliminate errors due to the finite size of the 
simulation cell. Details of the plane-wave cutoff 
and Appoint sampling will be given in the results 
section, where appropriate. 

The accuracy of DFT predictions depends 
mainly on the approximation used for the 
exchange-correlation functional. This is an im- 
portant issue in the present work, since the ad- 
sorption energy of hydrogen in Ca-GIC is fairly 
small - only a few hundred meV, and DFT er- 
rors are likely to be significant on this scale. To 
gauge the likely errors, we have performed the 
calculations with two rather different functionals: 
the LDA and the PW91 form of GGA. There is a 



general tendency for GGA functionals to be more 
accurate for binding energies, but their failure to 
produce any binding between the layers of pure 
graphite is a significant deficiency in the present 
context. On the other hand, LDA usually over 
binds, but its good prediction of the equilibrium 
interplanar spacing in graphite may be helpful 
here. The differences between LDA and PW91 
predictions can therefore be used as a rough indi- 
cation of the likely DFT errors. 

When searching for relaxed minimum energy 
structures, we generally use the automatic opti- 
misation algorithms of VASP to relax the ionic 
positions. The VASP algorithms can also be used 
to relax the size and shape of the simulation cell, 
but we found that it is often helpful to relax the 
cell parameters by hand, and we will note how 
this has been done at the appropriate places in 
the results section. 




Figure 1: The 6 starting positions of hydrogen in CaCg. 
Note positions (e) and (/) are both H2 molecules with (/) 
perpendicular to the page. 



3. Results 

3.1. Pure graphite 

As a check, we have repeated the calculations 
performed by earlier workers on pure graphite in 
the 2H, A — B — A structure that is observed ex- 
perimentally. We found that the LDA and PW91 
predictions of the in-plane C-C bond length are 
a=1.41 A and 1.43 A respectively, which are 
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very close to the experimental value of 1.42 A 
0, while LDA gives the value of c=3.38 A , 
in close agreement with the experimental value of 



n 



(1) 



3.36 A 14]. Our calculations on the A — A — A 
structure give the same values for a as in the 2H 
structure, and LDA gives c=3.60 A in this regis- 
tration. 

3.2. Calcium intercalated graphite (CaC n ) 

As noted above, we only need to study stage 
I Ca-GICs here. We shall present results for the 
different stackings a, a/3 and a/3y for CaCg but 
since the structural parameters and total energies 
turn out to be essentially identical, we examine 
only a stacking for other stoichiometries. The 
possible stoichiometries for CaC n are given by the 
formulae 

' 2 (k (k + 2) + 1) 
2(k(k — 1) + 1) 

with k = 1,2,3..., a positive integer. Here, we 
consider the cases k = 1,2,3, which correspond 
to n = 6,8,14. 

The unit cell was relaxed for all the Ca-GICs 
and the automatic algorithm in VASP used to 
relax the atomic positions within this equilib- 
rium cell. A plane- wave basis set energy cutoff 
of 800 eV and a /c-point grid with spacings always 
smaller than 0.04 A -1 were used for all calcula- 
tions. 

The two parameters that fully specify the struc- 
ture are the C-C in-plane bond length a and the 
inter-planar spacing c. Our calculated values of 
a and c are reported in table [lj which also in- 
cludes the available experimental data. We note 
that for CaCg, the differences between the differ- 
ent stackings are extremely small. The value of 
a is greater than that of graphite by only 1.5pm 
(1.0%) with PW91 and 2.5pm (1.8%) with LDA. 
The predicted values of c obtained with LDA and 
PW91 are fairly close to each other, and the latter 
is in excellent agreement with experiment. The 
value of c is greater than that of graphite by 
-34%. 

3.3. H 2 in CaC n 

Hydrogen was introduced into the CaC n struc- 
6, 8, 14, to test the effect of metal 



tures, with n 



density on the binding energy of H. By choos- 
ing only the the three stoichiometries mentioned 
we have limited the search for the equilibrium 
H position to be only within two inequivalent 
hexagons: one with the Ca atom in and one with- 
out. These are the only two distinct spaces that 
exist within CaC 6 , CaC 8 and CaCi 4 since every 
empty hexagon is equivalent. 

In order to sample these two spaces fully it was 
decided that hydrogen should be imported in both 
molecular form as H2 and atomically as 2H. Six 
distinct starting positions, shown in figure [TJ were 
chosen for hydrogen with positions (a) to (d) be- 
ing pairs of H atoms; and positions (e) and (/) 
being molecular H2, with all centres of mass lying 
perfectly mid-way between the graphite sheets. It 
should be noted that position (/) is perpendicular 
and position (e) is parallel to the graphite sheets. 
This results in six different atomic relaxations. 
Before these could be performed the hydrogenated 
CaC n unit cell was relaxed with all atomic posi- 
tions fixed in the starting position. Once the equi- 
librium starting structure was found an atomic 
relaxation was then started using the automatic 
VASP algorithm. A plane-wave basis set energy 
cutoff of 800 eV was used and a fc-point grid with 
spacings always smaller than 0.04 A -1 were used 
for all calculations. During relaxations (b) and 
(c) both relaxed back to position (a) showing the 
presence of a strong local minimum. 

The hydrogen bonding energies were calculated 
by subtracting the free hydrogen molecule energy 
and the relaxed CaC n energy from the hydro- 
genated CaC n , i.e. 

E (H 2 binding energy) = 

E (H 2 in CaC n ) - [E (H 2 ) + E (CaC n )] , 

thus giving a negative energy for favourable ab- 
sorption of H 2 . 

The results in table [2] show the absorption ener- 
gies for hydrogen in three different stoichiometries 
of CaC n in the most favourable position. In CaC 6 
position (a), two individual H atoms produced the 
highest binding energy. For CaCi4 position (/), 
molecular H 2 and perpendicular to the graphite 
sheets produced the highest binding energy. In 



Table 1: Calculated unit cell parameters for the CaC„ graphite intercalate and the comparison of CaC6 a/37-stacking 
with experiment 15 j. 





a(A) 


c(A) 




LDA 


PW91 


Expt. 


LDA 


PW91 


Expt. 


CaC 6 














a- stacking 


1.435 


1.445 


N/A 


4.34 


4.50 


N/A 


a/3-stacking 


1.435 


1.450 


N/A 


4.34 


4.47 


N/A 


a/?7-stacking 


1.434 


1.450 


1.444 


4.35 


4.50 


4.524 


CaCs a-stacking 


1.435 


1.445 


N/A 


4.31 


4.45 


N/A 


CaCi4 a-stacking 


1.424 


1.437 


N/A 


4.28 


4.45 


N/A 



CaCg the most favourable absorption was of an 
H2 molecule, with a negligible difference (approx- 
imately 3 meV) in binding energy between posi- 
tions ( e) (parallel to the graphite sheets) and po- 
sition (/) (perpendicular to the graphite sheets). 



Table 2: Calculated absorption energies for hydrogen in 
position (a) (two individual H atoms) for CaCg; and po- 
sition (/) (molecular H2 perpendicular to the graphite 
sheets) for CaCg and CaCi4. 





H 2 Binding 


Energies [eV] 




LDA 


PW91 


CaC 6 


-0.25 


-0.24 


CaC 8 


-0.08 


0.31 


CaCi4 


-0.01 


0.30 



It is interesting to compare this with hydrogen 
absorption in KC§ GIC For that material, the 
lowest energy configuration, by inelastic neutron 
scattering, is found to be one in which the H re- 
mains molecular and resides in a position near (/) 
jsj. Lovell et al. jij perform DFT calculations to 
find that the difference in H 2 binding energy be- 
tween the perpendicular (/) and the parallel (e) 
configuration in KCs is just a few meV, in quali- 
tative agreement with the predictions given here 
for CaCg. 

As shown in table |5] hydrogen was favourably 
absorbed with an LDA binding energy of 
—0.25 eV in CslCq. It is worth noting that when 
the relaxations were performed using the same 
plane wave cut off and k-point grid as Cobian 



and Iniguez 11] the results were in quantitative 
agreement with Cobian and Iniguez ll| . The hy- 
drogen was favourably absorbed with P W91 also, 
with a binding energy of —0.24 eV. The C-C bond 
length relaxed to a = 1.407 A, a modest —2.0% 
decrease upon pure CaC 6 . The inter-planar sep- 
aration changed to c = 5.24 A upon loading of 
hydrogen; a 20% increase, as expected from the 
formation of a Ca-H bond and thus the weakening 
of the Ca-C bond as noted by Cobian and Iniguez 
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3.4- Calcium hydride (CaH 2 ) crystal 

The equilibrium crystal parameters presented 
in this work are compared with the experimental 
results from Wu et al. 17J in table HI The cal- 



culated structure of calcium hydride also agrees 
well with the calculated structure of Li et al. [li| 
and Wolverton and Ozolins 19] . The agreement 



between calculated and experimental crystal pa- 
rameters is very good in PW91, and slightly worse 
with LDA. 



Table 3: Wyckoff 4 c coordinates within the CaH2 unit 
cell. Note that the x-, y- and z-coordinates are fractional 
coordinates of the lattice vectors. 



Position 


x-coord. 


y-coord. 


z-coord. 


1 


X 


0.25 


z 


2 


—x + 2 


0.75 


z+l 


3 


— X 


0.75 


—z 


4 


x + 2 


0.25 


-z+\ 
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Table 4: Calculated unit cell parameters and atomic positions for the CaH2 crystal and their comparison with experiment. 





Calculated K 




Experimental Ref. [17] 




(without zero point energies) 
LDA PW91 


9 K 


298 K 


t t 'i n 11 

Unit Cell 










Parameters 










a (A) 


5.72168 


5.90763 


5.92852(5) 


5.94753(6) 


6(A) 


3.46614 


3.57014 


3.57774(3) 


3.59326(3) 


c (A) 

XT / % 'iS 

V (A 6 ) 


6.56912 


6.77002 


6.78956(6) 


6.80185(7) 


130.280 


142.787 


144.011(1) 


145.362(2) 


Cal (4c) 










X 


0.23898 


0.23928 


0.2387(1) 


0.2397(2) 


z 


0.11180 


0.10988 


r\ "i "i r\c\ ( "i \ 

0.1102(1) 


0.1093(1) 


Til (C^^n \ 1 

ill ^CalC.J/ 
Dl (Exp.)(4c) 

X 


0.35479 


0.35586 


0.3558(1) 


0.3551(1) 


z 


0.42843 


0.42718 


0.4276(1) 


0.4268(1) 


H2 (Calc.)/ 
D2 (Exp.) (4c) 

X 


0.97416 


0.97482 


0.9750(1) 


0.9743(2) 


z 


0.68139 


0.67648 


0.6756(1) 


0.6759(1) 



To ascertain the stability of the loaded H 2 + 
CaC n system an energy comparison with CaH 2 
and graphite was performed. The stability was 
established by subtracting the energy of the H 2 
+ CaC n system from the separate energies of the 
calcium hydride crystal and graphite, i.e. 

AE = [E (CaH 2 ) + E (C n )] 

- E (H 2 in CaC n ) , (2) 

which gives AE to be negative if the H 2 + CaC n 
system is unstable against decomposition into the 
metal hydride and pure graphite. The calculated 
values of AE are reported in table El We see that 
both LDA and PW91 predict that all the loaded 
H 2 + CaC n structures are strongly unstable with 
respect to decomposition into graphite plus CaH 2 . 

4. Discussion and conclusions 

Our calculations completely confirm earlier pre- 
dictions that the intercalated graphitic material 
CaC 6 , does favourably absorb hydrogen, with an 



Table 5: Calculated stability energies, AE, against decom- 
position. 





H 2 + 
LDA 


AE 


[eV] for 
-> CaH 2 + C 6 
PW91 


AE for CaC 6 


-1.11 




-1.15 


AE for CaC 8 


-1.38 




-1.69 


AE for CaCi4 


-1.81 




-1.74 



absorption energy within the U.S. Department of 
Energy range of 0.2 — 0.8 eV. However, we have 
also found that the hydrogen-loaded CaC n mate- 
rial is strongly unstable with respect to decompo- 
sition into CaH 2 crystal plus graphite, the energy 
release associated with this decomposition being 
1.1 — 1.8 eV. This explains recent experiments 
12] , which showed that attempted hydrogenation 
of CaC6 prompted the complete deintercalation 
of the starting material to leave just CaH 2 and 
graphite. 



6 



In considering the reliability of our predictions, 
we must, of course, consider the likely size of er- 
rors with our DFT calculations. It is well known 
that DFT can easily suffer from errors of a few 
tenths of an eV in cohesive energies. We have 
seen in table [5] that the values of energy release 
for the decomposition of hydrogenated CaCg from 
LDA and PW91 differ by ~0.3 eV. A recent study 
of the energetics of MgH2 by Pozzo and Alfe 20 
showed that DFT errors in the cohesive energy 
can be as large as 0.5 eV, but that LDA and 
PW91 predictions bracket the true value. How- 
ever even with errors of this size, our calculations 
leave no doubt that the decomposition is strongly 
exothermic, with a decomposition energy of at 
least 1.0 eV. 

It is worth noting that thermodynamic instabil- 
ity with respect to decomposition does not neces- 
sarily tell us that decomposition will occur spon- 
taneously at room temperature, because the de- 
composition might conceivably be kinetically hin- 
dered. However, the experiments 12| show that 
at room temperature and within one hour of hy- 
drogen exposure the sample lost all evidence of 
the CaC6 structure, and within five hours showed 
only the CaH 2 and graphite structures. These 
time-scales suggest that the rate-limiting step 
probably has an activation energy of ~0.9 eV. The 
computational search for possible decomposition 
pathways would be a major undertaking, and is 
well beyond the scope of the present work. 

The results we have found raise the interest- 
ing question of what determines whether a hy- 
drogenated metal GIC will be unstable or not 
with respect to decomposition. We know that 
hydrogenated Li-GIC and Ca-GIC are both un- 
stable. However, K-GIC appears to be stable, 
and a number of experiments have been reported 

on the hydrogenated material QBElil. we 

suggest that the dominant factor may be the co- 
hesive energies of the metal hydride, with those 
having large cohesive energies being unstable. 
This would probably fit the facts, because the 
larger lattice parameter of KH compared with LiH 
should give it a smaller cohesive energy. The dou- 
ble charge of Ca in CaH 2 should also make its co- 
hesive energy larger than that of KH. This is an 



interesting question for future work. 
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